Streamwise-traveling waves 

of spanwise wall velocity 
for turbulent drag reduction 



MAURIZIO QUADRIO , 
PIERRE RICC0 2 and CLAUDIO VIOTTI 3 

I Politecnico di Milano, 20156 Milano, Italy 
2 King's College London, WC2R 2LS London, United Kingdom 
3 The University of North Carolina, 27516 Chapel Hill, NC, USA. 

(Received 8 January 2009) 

Waves of spanwise velocity imposed at the walls of a plane turbulent channel flow 
are studied by Direct Numerical Simulations. We consider sinusoidal waves of spanwise 
velocity which vary in time and are modulated in space along the streamwise direction. 
The phase speed may be null, positive or negative, so that the waves may be either 
stationary or traveling forward or backward in the direction of the mean flow. Such a 
forcing includes as particular cases two known techniques for reducing friction drag: the 
oscillating wall technique (a traveling wave with infinite phase speed) and the recently 
proposed steady distribution of spanwise velocity (a wave with zero phase speed) . 

The traveling waves alter the friction drag significantly. Waves which slowly travel for- 
ward produce a large reduction of drag, that can rclaminarizc the flow at low values of the 
Reynolds number. Faster waves yield a totally different outcome, i.e. drag increase. Even 
faster waves produce a drag reduction effect again. Backward-traveling waves instead 
lead to drag reduction at any speed. 

The traveling waves, when they reduce drag, operate in similar fashion to the oscillating 
wall, with an improved energetic efficiency. Drag increase is observed when the waves 
travel at a speed comparable with that of the convecting near-wall turbulence structures. 
A diagram illustrating the different flow behaviors is presented. 



1. Introduction 

Reducing the friction drag in turbulent wall flows has seen a growing interest in recent 
years, owing to the awareness for high energy consumption and to the need of decreasing 
pollutant emission s into the atmosphere. Drag reduction techniques are usually grouped 
( Gad-el Hak 2000| ) into active (requiring external energy to be fed into the system) 



and passive techniques. Another useful classification is between closed-loop strategies 
(requiring a control law with feedback) and simpler open-loop techniques: the review by 



Kiml ((2003) discusses the former group, whereas the present work is concerned with the 
latter. 

One interesting open-loop technique, which is conceptually simple and with the poten- 
tial for a positive energetic budget, is the cyclic spanw ise movement of the wall. In this 
active technique, first introduced by I Jung et al. ( 1992f ). the wall moves sinusoidally with 



period T according to 

w w (t) = Asm (ujt) , (1.1) 



2 M. Quadrio, P. Ricco & C. Viotti 

where w w is the spanwise (z) component of the velocity vector at the wall, t is time, A is 
the oscillation amplitude and lo = 2tt/T is the oscillatio n frequency. The main features 
of the oscillating- wall technique ( Ricco fc Quadrioll2008f ) are the existence of an optimal 



frequency Lo opt for drag reduction at fixed A, and a significant maximum drag reduction, 
which can be as high as 45% when A is comparable with the flow centreline velocity. At 
lower A, the energy balance between the power saved thanks to the oscillation and the 
power spent to activate the wall motion against the viscous stresses may be positive at 

about 7% when lo = L0 njlt . 

In a recent work, Quadrio et all ( 20071 ) (see also Viotti et al. 120081 ) have addressed 



one of the drawbacks of the oscillating-wall technique, i.e. its unsteady nature. They 
converted the time-dependent forcing (|l.lj) into its stationary counterpart of a spatial 
sinusoidal distribution of spanwise velocity over a wavelength A^: 

w w (x) — Asm (k x x) , (1-2) 

where x is the streamwise coordinate and k x = 2n/X x is the streamwise wavenumber. The 
time dependence of Eq. (|1.1[) is converted into a relatively unsteady interaction between 
the steady wall motion and the convecting near-wall turbulence. Such conversion is only 
apparently similar to the Taylor's frozen-turbulence hypothesis, for which the velocity 
scale is the mean streamwise velocity. Indeed, the mean velocity is zero at the wall , 
while the turbulence keeps a strongly convective character ( Kreplin fc Eckelmannl 19791 : 



Kim fc HussainllT993l ). with elongated space-time correlations ( Quadrio fe Luchinill2003h . 



The velocity scale for space-time conversion is therefore the convection velocity 11 of the 
turbulent fluctuations, which almost coincides with the mean velocity along with the 
wall-normal span of the channel, except for a thin near-wall layer, say y + < 10, where it 
is constant at 10. (Hereinafter, the + superscript indicates scaling b y the friction 

veloci ty u T of the reference flow and the kinematic viscosity v of the fluid.) IViotti et 



( 20081) ' have shown that the optimal frequency 0J O pt for (|1.1[) translates into an optimal 
wavenumber K x , op t for (| 1 . 2[) through the relation 

K x ,opt = • (1-3) 

The present paper addresses the extension of (|1.1|) and (II. 2|) to the space-time case 

w w (x, t) = Asm (k x x — tot) , (1-4) 

where n x and lo may be both different from zero at the same time. In p.4[) . the wave 
of spanwise velocity moves (backward or forward) in the streamwise direction, and its 
phase speed is 

lo 

c=— . (1.5) 

k x 

The temporal (|l.ip and steady (|1.2[) oscillations correspond respectively to a wave trav- 
eling at infinite speed (k x — 0), and to a stationary wave (lo = 0). The effect of the wall 
motion is already known on the axes of the lo — k x plane, while it is unexplored in the 
rest of the plane. The main question addressed in this paper is therefore: what is the 
response of the turbulent flow in the whole lo — k x plane? 

Before presenting our resu lts, it may be u s eful t o mention the drag reduction technique 
bv lDu fc Karniadakid (|2000h and lDu et al\ (|2002h . They studied the following spanwise- 



oriented volume forcing: 

f z {z, t) = F e-y/ A sin [n z z - Lot) , (1.6) 



where A is the distance up to which the forcing diffuses from the wall and k z is the 
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spanwise wavenumber. In a turbulent channel flow at Re T — 150, they obtained about 
30% drag reduction and observed a substantial disruption of the near-wall flow struc- 
tures. The volume force in Eq. (|1.6p consists of spanimse-traveling waves acting along the 
spanwise direction, while here we consider a wall-based motion of spanwise velocity waves 
traveling along the streamwise direction. 

In the present paper, the effects of imposing condition (|1.4p are investigated by Direct 
Numerical Sim ulations (DNS). W hile purely spatial and purely temporal oscillations have 
been shown bv lViotti et al. (2008) to be largely analogous to each other, the flow response 
to the traveling waves is found to be interestingly complex when examined as a function 
of u> and k x . Both drag reduction (DR) and drag increase (DI) may occur. The maximum 
DR resides in a non-obvious region of the lo — k x plane, which largely corresponds to 
where the global energy balance is positive, whereas the maximum DI occurs where the 
phase speed c is comparable with the convection velocity IA W . 

The structure of the paper is as follows. In <J3J the numerical procedures are outlined. 
Section Sj3]presents the drag modification and the global energy budget in the u>—k x plane, 
and discusses the effect of forcing amplitude and Reynolds number on the maximum DR. 
Selected flow statistics arc also described. In fjU the analogy between the oscillating wall 
and the traveling waves is discussed, and a diagram offers a comprehensive view on the 
modified turbulent flow. A brief summary is given in Sj5] 



2. Method 

A large number of DNS has been carried out to examine the response of a turbulent 
channel flow to the wall forcing (|1.4[) . The traveling waves are applied at both walls 
and move in phase. From a numerical perspective, Eq. (|I ,4p is a non-homogeneous, time- 
dependent wall boundary condition for the spanwise component of the velocity vector. 
The size of the computational domain, the requirements in terms of the spatio-temporal 
resolution and the length of the integration interval render these DNS very demanding. 
The efficiency of the parallel DNS code and the availability of appropriate computational 
resources are key aspects of this research effort. 

The code and the architecture of the computing system have been described bv lLuchini fc Quadrio 
(l2006h . The code is a mixed-discretization parallel solver of the incompressible Navier- 
Stokes equations, based on Fourier expansions in the homogeneous directions and high- 
order, explicit compact finite-difference schemes in the wall-normal direction. Its primary 
design target is to achieve high performance while using off-the-shelf hardware. Most of 
the calculations have been run on a computing system available in dedicated mode at the 
Universita di Salerno. The system is composed of 268 AMD Opteron CPUs, grouped in 
134 computing nodes, and connected with an ad hoc topology through Gigabit Ethernet 
cards. Each node is equipped with its own hard disk and RAM. The system possesses 0.28 
TB of RAM and 40 TB of storage space, with a peak computing power of 2.6 TFlop/s. 

The computations are integrated forward in time, starting from the same initial con- 
dition of a fully developed channel flow without waves. The flow rate is constant and the 
bulk velocity is £/& = 2/3£/p, where Up is the centreline velocity of a laminar Poiseuille 
flow with the same flow rate. About 250 simulations have been run with a wave amplitude 
A = 0.5Up (which corresponds to A + = 12) at a Reynolds number Re = Uph/v = 4760 
(which corresponds to Re T — hu T /v = 200 based on u T and the channel half width h). A 
few cases have also been run to explore how the maximum DR is sensitive to changes in 
A or Re: A/Up has been varied in the range 0.085 — 1.26, and two values of Re = 2175 
and Re = 10500 (corresponding to Re T = 100 and Re T = 400, respectively) have been 
considered. 



4 



M. Quadrio, P. Ricco & C. Viotti 




Figure 1. Schematic of the system for turbulent channel flow with wall traveling waves. 
is the streamwise wavelength and c is the phase speed of the waves (traveling forward in this 
sketch). L x , L y and L z are the dimensions of the computational domain in the streamwise, 
wall-normal and spanwise directions, respectively. 



The computational domain, a schematic of which is presented in figure[I] has a stream- 
wise length of L x = Qirh and a spanwise width of L z = inh. At Re = 4760, a Fourier 
expansion of 320 x 320 modes is used in the streamwise and spanwise directions, and 
160 points discretize the wall-normal direction, for a total of 3.3 • 10 7 spatial degrees of 
freedom. The streamwise and spanwise resolutions are Aa; + = 11.8 and Az + = 5.9, and 
the wall-normal resolution varies from At/^ = 0.7 at the wall to Ay+ aa . = 4.1 at the 
channel centreline. The time step is At + = 0.16, and the total integration time is 1000 
h/Up (i.e. 8400 viscous time units). A few cases (first row of data at K x h = 0.167 in 
figure [2]) have L x = 12nh, and use 640 streamwise Fourier modes to keep the resolution 
unchanged. A few other cases (with extremely large oscillation periods) have an averag- 
ing time of 2000 h/Up. For the simulations at different Re, the size of the computational 
domain has been kept constant in outer units, and the spatial resolution has been kept 
constant in viscous units. The effect of the discretization para meters on the fricti o n dra g 
has been assessed by preliminary checks. In agreement with Quadrio fc Riccol ( 2004 ). 
the accuracy in measuring the friction coefficient is better than 1%, and can be further 
improved by extending the total averaging time. 

The main objective is to investigate how the forcing (|1.4p alters the friction coefficient, 
defined as 

Cf = W b ' 

where t w is the mean wall-shear stress and p is the density of the fluid. Changes of 
Cf are quantified in terms of percentage of C/ ; o, the friction coefficient of the reference 
turbulent flow, and coincide with the percentage saving of the energy required to drive 
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the fluid along the x direction at fixed flow rate. For the reference case at Re = 4760, 
Cf.p = 7.94 • 10~ 3 , which is in very good agreement w ith the correlation reported by 
Popd (|2000h at page 279. As in lQuadrio k Riccol (|2004h . C f is evaluated by space- and 



time-averaging, after discarding the initial transient during which the flow adapts to the 
new regime. The length of the transient is estimated by direct observation of the time 
history of the space-averaged friction, and is typically of the order of 100-200/i/C/ p . Unless 
otherwise indicated, quantities reported hereinafter are scaled by h and Up. 



3. Results 

3.1. Effect on friction drag 

Figure^] describes the effects of the traveling waves in terms of the percentage change in 
friction drag as a function of u and k x for A — 0.5 and Re = 4760. Only the upper half of 
the u) — k x plane is shown, since (jl.4p yields symmetric results upon exchanging the pair 
u>, k x with — u>, ~n x . (This has been explicitly verified for a few points.) The phase speed 
c is interpreted graphically as the inverse of the slope of straight lines passing through 
the origin. In the first quadrant, the waves move forward in the direction of the mean 
flow (c > 0), while the second quadrant corresponds to backward-traveling waves (c < 0). 

Each point in figure O is the result of one simulation. The contours assist in visualizing 
the general behaviour of the changes in drag; they are computed after a linear interpo- 
lation of the irregularly scattered points on a regular mesh with grid sizes of Aw = 0.05 
and Ak x = 0.1. The presence of the numerical values, rounded to the nearest integer, 
allows verifying that the contours do not misrepresent the data. Several points with large 
uj and/or k x do not appear, but they are however useful to obtain the correct trend of 
the contours near the boundaries. 

On the horizont al uj axis. Eg. (11.41) reduc es to the oscillating- wall case (jl.ip : our data 



agree with those bv lQuadrio fc Riccol ((20041), who used the same Re and slightly different 



discretization parameters. The symmetry mentioned above requires points at uj and — uj 
to have the same friction drag: this is enforced by plotting every data point twice, at 
+ui and — uj. On the vertical k x a xis. Ea.(|1.4l) reduce s to the steady wave case (|1.2p . 
Some of these points are taken from lViotti et all ( 20081 ). who employed the same Re and 



discretization parameters. As discussed there, DR behaves very similarly on the two axes, 
once frequency is converted into wavenumber through Eq. (|1.3[) . reaching its maximum 
at ui = uj opt k, 0.5, and at k x = n x , pt ~ lj with uj opt and n Xt0pt related by U w = 0.5. 

Off the axes, the flow response reveals an unexpectedly rich behaviour. The contours 
show that local maxima and local minima form two elongated narrow regions. In the red 
region, whose crest identifies a curve which does not cross the origin, the flow presents 
an intense drag reduction. The maximum is about 48%, which is higher than 34% for 
the oscillating wall at ui opt and 45% for the steady wave at n x _ opt . A wide range of k x 
exists where DR is 47%-48% for a non-zero phase speed. The crest crosses the vertical 
axis at k x = K x ^ opt and intersects the uj axis at uj — —ui opt . For k x > K Xt0p t, this region 
largely corresponds to slow forward-traveling waves (c « 0.15), while for k x < K x ^ opt it 
relates to backward-traveling waves. In the blue region, the drag increases significantly, 
up to 23%. The maximum increase corresponds to a straight line crossing the origin, and 
is therefore associated with a constant phase speed, c = 0.5. The narrow, cone-shaped 
DI region is confined between the two thick zero-DR lines, corresponding to c 0.35 and 
c « 0.6. The marked asymmetry between the first and second quadrant vanishes at high 
\uj\: the contour lines tend to become vertical, suggesting that in this limit the effect of 
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Figure 3. Friction drag reduction (percentage) as a function of T = 2-k/lu for A = 0.5, 
Re — 4760 and k x — 0. The leftmost points are at T = 1 and T = 1.4, and points with 
2 < T < 50 are out of scale. 



the streamwise modulation becomes negligible and to alone dictates the drag reduction 
level. 

Careful observation of the region near the origin of the u> — n x plane reveals a finite 
DI for the oscillating wall at very small fr equencie s. This is shown in figure EJ where 
T = 2it/lu is used as independent variable. Nikitin ( 200fj| ) was the first to surmise the 
existence of what he cal led an infrared DI, based on an under-resolved DNS at T > 120. 
Quadrio fc Riccol ( 2004 ) tested the oscillating- wall technique in a plane channel up to 
T k 90 and were unable to confirm the infrared DI (figure [3] indeed confirms that no 
drag increase is observed for T < 100). Owing to the computational cost, we did not 
explore this issue further. 

3.2. Power budget 

The efficiency of the traveling waves in affecting the friction drag is studied by computing 
the amount of external energy spent to enforce the control action. The quantification of 
this energy obviously implies the knowledge of the actuation device. However, the power 
that an ideal device dissipates against the viscous stresses provides an upper bound for 
the performance of the control strategy. We express this power as percentage of the power 
spent to drive the fluid along the streamwise direction in the fixed-wall configuration: 



P sp (%) = 



100 



1 



(dU/dy\ )U b L X L Z T 



dw 
' dy 



dtdzdx. 



where dU/dy\o is the gradient of the mean streamwise velocity at the wall in the reference 
case. P sp (%) is reported in figureS]as a function of to and k x . Much less energy is required 
to impose a standing wave than to generate uniform oscillations, if cases related through 
Eq. (|1.3p are compared. However, neither the oscillating wall nor the standing waves are 
associated with the pair which offers the minimum energetic expense of P sp (%) = 23.5, 
i.e. to w 0.15, k x ~ 1. 
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Figure 4. Map of the power P sp (%) spent for the traveling waves in the ui — k x plane for A = 0.5 
and Re = 4760. Contours are spaced by 10% intervals and the minimum is 30%. Here and in 
the following figure, the percentage is relative to the reference friction power required to drive 
the flow along the streamwise direction. The numbers indicate P ap (%) at measured points. 



Even more important is to quantify the net energy saving, defined as the difference 
between the percentage power saved to drive the fluid along the streamwise direction 
(which coincides with the drag reduction at constant Ub) and the percentage power 
spent defined above, namely 

Pnet(%) = DR(%) - P sp (%). 

For several acti ve open-loop techniques , P ne t(%) is negative. Some methods, like the 
oscillating wall ( Baron fc Quadriol 19961 ). may offer a positive net budget of the order 



of a few percent. The traveling waves, on the other hand, guarantee much larger ben- 
efits because of the remarkable occurrence that the region of minimum P sp (%) largely 
coincides with the region of maximum DR(%). Figure [5] shows that P ne t{%) is positive 
mainly for slow forward-traveling waves, and that a 18% of maximum net is measured 
for lu w 0.15, k x ~ 1 (which is also the optimum pair for minimum P sp ). This benefit 
is noteworthy when compared with the negative net balance of the oscillating wall at 
this value of A, and the small net gain of 5% obtained by the steady waves. It must be 
recalled that we have discussed a simple estimate, which cannot account for the nature 
of the real actuation system. However, such a high positive net power budget might well 
accommodate a real system with efficiency less than unity. 
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Figure 5. Map of net power saving P ne t(%) in the ui — k x plane for A — 0.5 and Re = 4760. 
Contours are spaced by 10% intervals. Locus of zero P n et(%) is indicated by thick curve, solid 
lines denote positive balance and negative values are represented by dashed lines. The numbers 
indicate P„ e t(%) at measured points. 



3.3. Effect of A and Re on maximum drag reduction 

The sensitiveness of the maximum DR to the parameters A and Re, kept fixed in the 
main parametric study, is addressed through a few additional simulations. Figure [5] (left) 
shows that DR increases monotonically with A and reaches nearly 60% when A is slightly 
larger than the flow centreline velocity. In the plot on the right, P ne t grows significantly 
at lower amplitude, which is expected as the power spent to impose the traveling waves 
decreases quadratically with A, and is larger than 26% for A = 0.25. Both results are con- 
sis tent with (but quantitat ively much better than) the oscillating-wall scenario, studied 
bv lQuadrio fc Riccol (|2004) . 

The maximum DR at Re = 4760 depends weakly on the Reynolds number. We have 
observed full relaminarization of the flow when Re = 2175 (Re T = 100), and a slight 
decrease from 48% to 42% as Re increases to 10500 (Re T = 400). T his finding is con- 
sistent with experiment al ( Choi fc Graham 1998; Ricco_& Wull2004h and DNS results 



( Ricco fc Quadrioll2008l ) available for the oscillating-wall case. Following these works, we 



have assumed that the values of lo and k x that guarantee the maximum DR scale in wall 
units of the reference flow. 

3.4. Turbulence statistics 

In this section, we study the turbulence statistics to understand how the flow is modified 
by the traveling waves. We first take a preliminary look at instantaneous snapshots of the 
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Figure 6. Percentage drag reduction (left) and net power saving (right) as functions of forcing 
amplitude A. The forcing conditions are u — 0.16, k x = 1.66, which produce the maximum drag 
reduction at A = 0.5 and Re — 4760. 

flow field in figure [7] Two different cases are compared to a snapshot of the reference flow, 
shown in the top plot: one with u) = 0.16, k x — 1.66 and the other at lu — 0.88, n x = 1.66. 
These two cases do not lie far from each other in the lo — k x plane, but yield totally 
different outcomes, i.e. a large DR of 48% for the former (middle plot) and a large DI 
of 14% for the latter (bottom plot). Both cases refer to forward-traveling waves, with 
c = 0.1 and c = 0.53, respectively. 

Isosurfaces for negative values of the quantity q = sgn(u)-\At 2 + w 2 are shown in figure 
[7] to visualize the low-speed streaks of the near- wall flow even in the present case where 
the local streamwise direction may deviate from the x axis. The value q = —0.168 (q + = 
—4) is chosen for the reference flow, that indeed shows the usual pattern of near-wall 
elongated, streamwise- aligned regions. In the plots for the DR and DI cases, the value of q 
is changed so as to always yield q + — —4 when computed with the actual friction velocity 
of the flow. This compensates for the largely different values of Re T : had the nominal 
value of u T been used, the streaky structures would have disappeared in the DR case. 
When DR takes place, what remains after accounting for the reduced Re T are elongated 
structures, aligned with the x direction, of perhaps larger size and reduced streamwise 
meandering. In sharp contrast, in the DI case the same structures present an evident 
streamwise modulation, with a streamwise wavelength corresponding to the wavelength 
of the wall wave, which is a marked structural difference with the reference flow. In this 
case the most evident action of the forcing is to cyclically change the direction of the flow, 
by tilting the low-speed streaks in accordance with the local spanwise velocity component 
at the wall. 

In figure[51 the profiles of the mean streamwise velocity component, averaged over time, 
x and z, are presented as functions of the wall-normal coordinate. They are indicated by 
the solid line (reference case), filled (DR) and empty (DI) circles. The effect of the wall 
motion diffuses all along the extent of the channel and the large near-wall change of the 
velocity gradient is evident when DR occurs. 
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Figure 7. Snapshots of instantaneous flow fields for the reference flow (top), a case with DR 
(middle, n x = 1.66 and w = 0.16) and a case with DI (bottom, n x = 1.66 and uj — 0.88). 
Isosurfaces for the quantity q = sgn(u)V'« 2 + w 2 at the value q + = —4. The value q + = —4 is 
computed with the actual friction velocity of the drag-modified flow. The flow is from left to 
right and only the bottom half of the channel is shown. For both cases A = 0.5 and Re — 4760. 
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Figure 8. Profiles of mean streamwise velocity U, after averaging over time, x and z. The solid 
line is the reference case, the empty circles are the DI case (uj = 0.88, K x = 1.66), and the filled 
circles are the DR case (uj = 0.16, k x = 1.66). For both cases A — 0.5 and Re = 4760. 




Figure 9. Contour lines of the (w) Zt t field. Contours are (—0.45,0.1,0.45), dashed lines cor- 
respond to negative values. In this figure and in the following ones, the graphs are for one 
streamwise wavelength. The left graph refers to the DI case (uj = 0.88, k x = 1.66), while the 
right one refers to the DR case (uj — 0.16, k x = 1.66). For both cases A = 0.5 and Re = 4760. 



The mean spanwise velocity and the components of the Reynolds stress tensor are 
now considered. The averaging is over time and z, at different oscillation phases over 
the convective variable £ = x — ct, which is the proper phase alignment with the wall 
motion. (Hereinafter, this averaging operator is indicated by (-)z.t-) In figures [51 and [TU1 
the quantities are represented by contour plots as functions of £ and y. The streamwise 
direction is from left to right and for each case the left plot refers to the DI case, while 
the right contour indicates the DR case. The computational domain spans a streamwise 
distance of five wavelengths, but only quantities corresponding to one wavelength are 
shown. 

Figure [H displays the mean spanwise velocity, i.e. (w) z j- Although the effect on the 
friction drag is opposite in the two cases, the appearance of the two fields is rather 
similar. A phase lag of spanwise velocity along the vertical direction is observed: at a 
non-zero wall-normal location, the value of the wall velocity is encountered at a different 
£ position. Similarly to the case of uniform oscillations, the convecting flow over the 
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traveling waves is exposed to an unsteady transversal boundary layer. The thickness of 
the two ^-modulated boundary layers in figure [9] is comparable. 

The six components of the Reynolds stress tensor, i.e. (uiUj) z ,t, are shown in figure 
[TU1 As expected, the intensity of all the components consistently increases in the DI 
case (when compared to the reference flow, where these quantities are a function of y 
alone) and is significantly attenuated in the DR case. All the quantities in the DI case 
except (vv) Zt t show a marked modulation along £, consistently with the instantaneous 
snapshot of figure [7] (bottom). They can be qualitatively connected to the statistics of 
the reference flow, if the locally changing flow direction is accounted for. The near-wall 
peaks of (uu) Zjt , for example, are related to similar peaks of (ww) z j, at different phases. 

The variation along £ is less intense in the DR case. The reduction of the levels of 
{uv) z ,t, pertaining to the mean flow effects on the wall friction, is particularly noteworthy. 
The peak of the (uu) z ,t contours is closer to the wall in the DI case than in the DR case, a 
well-known phenomenon in the oscilla ting- wall flow, for which the viscous and the buffer 
layer thicken when friction decreases ( Choi et al. 1998). For the (vw) Zl t component, the 



level of fluctuations is so small in the DR case that no contour is shown on the chosen 
contour scale. Note moreover that (uw) z t and (vw} z t are zero in the reference flow. 



4. Discussion 

The mechanism by which the traveling waves affect the turbulent friction is described 
in this section. We surmise that this mechanism includes the one of the oscillating wall, 
complemented by an additional effect responsible for drag increase and due to the phase 
speed attaining a finite value. 

This additional effect is linked with th e fact that the turbulen t flow at the wall possesses 
the convective velocity scale U w — 0.5 (jKim fc Hussainlll993f h Indeed, when the phase 
speed of the waves approaches U w , a large DI occurs. This may be interpreted by a phase 
lock between the near-wall flow structures and the wall waves, which is strongest along 
the line of maximum DI in figure [21 i.e. when c = U w . Since U w is a single velocity value 
computed from flow statistics, it should be viewed as representative of the (different) 
convection speeds of the near- wall structures, which can be smaller or larger than U w . 
Consistently with this interpretation, the triangular region where DI occurs is identified 
by a range of phase speeds. When DI occurs, the structures are cyclically tilted by the 
wall forcing, which is stationary when observed in a reference frame convecting with the 
structures. The resultant streamwise modulation of the structures in the DI case is well 
illustrated in figure [7J 

Away from the triangular region, the mechanism by which the traveling waves influ- 
ence the friction drag becomes similar to that of the oscillating wa ll. In this latter case 
of un iform motion, the generally accepted qualitative explanation ([Karniadakis fc Choi 
2003h is that, when the wall oscillates in time with a period comparable with an opti- 



mal T opt = 2tt /ijj op t, the transversal Stokes layer effectively interacts with the near- wall 
turbulent structures, weakens the near- wall viscous cycle and produces a large DR. The 
optimal period is related to a time scale re presenting the typical life time of the most 
statistically significant turbulent structures (jQuadrio fc Luchinill2003l ). and to an opti- 
mal thickness of the Stokes layer. At periods smaller than T opt , the Stokes layer is thin 
and the interaction with the turbulent flow is limited to the very near-wall region. DR 
increases with T at these small periods. When T is larger than T opt , the Stokes layer 
becomes too thick and the interaction with the turbulent flow is not as effective: the DR 
effect rapidly deteriorates for increasing T. 

This general picture can be easily adapted to the traveling waves. A first occurrence 
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Figure 10. Components of the Reynolds stress tensor. From top to bottom: {uu) z ,t, {vv) z ,t, 
{ww) z ,t, (uv) z ,t, (uw) z ,t and {vw) z ,t- Contours are (0,0.0005, 0.01) for the diagonal components, 
and (—0.0042, 0.0004, 0.0042) for the off-diagonal components; dashed lines correspond to nega- 
tive levels. The graph on the left refers to the DI case (ui = 0.88, n x = 1.66), while the graph on 
the right refers to the DR case (w = 0.16, k x = 1.66). For both cases A = 0.5 and Re — 4760. 
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that points toward such a generalization is that the thickness of the transversal boundary 
layer assumes very similar values when the oscillating wall, the standing wave and the 
traveling wave guarantee the maximum DR. If the thickness is defined as the wall- normal 
location where the maximum (over all phases) of the spanwise velocity reduces to e _1 A, 
the optimal thickness corresponding to maximum DR is 3.15 • 1CT 2 for the oscillating 
wall, 3.6 • 10 -2 for the steady waves, and 3.4 • 10~ 2 for the traveling waves. 

It can therefore be conjectured that in all the three cases the distance of viscous 
diffusion from the wall toward the core of the turbulent flow is important in determining 
the turbulence suppression. The key difference with the oscillating wall is that the period 
to be compared with T opt is not the oscillation period T of the traveling wave, but an 
equivalent period T which is related to the convective nature of the near-wall turbulent 
fluctuations, and to both w and n x . We define T as the time required by the near- 
wall turbulent structures to cover a distance of one wavelength along the streamwise 
direction at the relative velocity IA W — c. The simple formula for T is 



T is thus the time interval between two consecutive instants at which a near-wall structure 
convecting at velocity U w passes over the waves at the same phase. 

Analogously to the oscillating wall, two cases can be distinguished. When |T| < T opt , 
or, equivalently, when 



the forcing due to the traveling wave yields a thin transversal boundary layer, the bulk 
flow remains aligned along the streamwise direction, and DR at constant k x increases 
with |T|. When \T\ > T opt , i.e. when 



DR at constant n x instead decreases with |T|. 

In the ui — k x plane, the points satisfying (|4.2[) identify the exterior of the oblique strip 
bounded by two parallel straight lines with slope k x /uj — (Uw)^ 1 (and thus parallel to 
the line of maximum DI) passing through the points to — ±oj opt , k x — 0. The condition 
T = T opt identifies the maxima at to = 0, k x = n Xt0p t = 1 for the standing waves, and 
at fc x — 0, u) = Lo op t = —0.5 for the oscillating wall. For backward traveling waves, the 
analogy with the oscillating- wall regime holds. Moving along lines of constant k x , DR 
increases with T, except, as expected, inside the strip where DR decreases with T. The 
red ridge of maximum DR in figure [2] and the border of the strip T = T opt show good 
agreement. 

For forward traveling waves, the behaviour of DR is determined by the mutual action 
of the two mechanisms outlined above. The ridge of maximum DR now lies entirely at 
the exterior of the strip and not at the border of it. This is because the DI region widens 
at large k x and dominates the oscillating-wall condition to dictate the changes in drag. 
Points laying outside the strip and outside the influence of the triangular region show 
increasing DR at fixed k x for increasing \T\. At small n x , points at low frequency are 
inside the strip where DR drops as \u)\ decreases, ultimately falling into the infrared DI 
region. 

The sketch in Figure [IT] summarizes the effects of the traveling waves on the turbulent 
flow. The thick solid curve denotes the locus of points of maximum DR at given u>. The 
friction drag decreases everywhere in the plane, except in the cone-shaped area, where DI 
occurs. The dash-dotted line indicates the points of maximum DI at given k x , occurring 



\14 W K X — Uj\ > L0 op t, 



(4.2) 



\14 W K X — Ul\ < L0 op t, 



(4.3) 
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Curve of maximum 



Maximum DR 
(traveling waves) 



Maximum DR 
(standing waves) 



Line of maximum DI 




Maximum DR 
(oscillating wall) 



Figure 11. Qualitative sketch of the map of DR due to the traveling waves in the uj — k x 
plane. DR occurs in the entire plane, except in the dashed area where DI takes place. Thick 
lines indicate loci of maximum DR (solid) and maximum DI (dash-dotted) at given n x . The 
region confined between the two oblique dashed lines is where the wall forcing is slower than 
the optimal timescale, as observed in the moving reference frame (see text). Symbols indicate 
points of maximum DR on the axes (closed circles) and on the entire plane (open circle). 



when c = U w , or, equivalently, when T — > oo. The dashed lines mark the strip defined 
by|T|=T opt . 

4.1. Open questions 

Several interesting questions are left unanswered by the discussion above. A non-exhaustive 
list includes: Can the details of the DR mechanism be better understood? What are the 
real possibilities of this technique in terms of net power savings? Do particular wave- 
forms exist, possibly different from sinusoids, that allow for further optimizations? There 
is clearly need for a further understanding of the mechanisms by which the wall friction 
is modified by the traveling waves. 

An additional, important issue is how to exploit the present results for practical pur- 
poses. This is unfortunately not immediate, but a sound understanding of the com 



plex fl ow physics is definitely an essential preliminary step. As discussed by IViotti etal 



(2008), the discovery of a preferred longitudinal length scale selected by the standing 



waves could have a significant impact on the design of optimal wall roughness for drag 
reducti on. Recent works have addressed the application of similar traveling- waves tech- 



niques. lZhao_ei_aL| (|2004[ ) have pointed out that an array of smart materials and MEMS 



would render a smart flexible wall the most promising candidate for effective turbu- 
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lence control. Itoh et al. ( 20061 ) have succeeded in studying experimentally the spanwise- 
traveling wave concept in a turbulent boundary layer through a flexible sheet. 

One further point of interest is the prediction of the drag modification (without 
resorting to direct measurements either via full DNS or laboratory experiments). In 
the oscillating-wall case, the space-averaged velocity profile follows cl osely the laminar 
solution, which has le d to a scaling parameter f o r dra g reduction (jChoi et all 12002 ; 
Quadrio fc Riccol 120041 ). Analogously. IViotti et all §008) have found that the laminar 
layer induced by the standing waves agrees well with the space-averaged turbulent profile. 
Future work should therefore be directed at extending these results to the traveling-wave 
case to estimate the modification of drag. 



5. Summary 

This work has described the response of a turbulent channel flow to sinusoidal waves 
of spanwise velocity applied at the wall and traveling along the streamwise direction. 
The effects of the wavenumber k Xi temporal frequency to and - to a lesser extent - 
forcing amplitude A and Reynolds number Re have been explored via Direct Numerical 
Simulations. Such a demanding parametric study, that required centuries of CPU time, 
would have been impossible, had our special computing system not been available. 

In the lu—k x plane, a region where drag reduction (DR) is very large has been identified. 
This region pertains to forward-traveling waves with small phase speed, but extends in 
the region of backward-traveling waves at small k x . A maximum DR of 48% is measured 
for Re = 4760 and A — 0.5Up; DR above 45% is observed over an extremely wide range 
of wavelengths, provided the frequency is adjusted to maintain the correct phase speed. 
DR weakly decreases with the Reynolds number, going from full relaminarization of the 
flow at Re = 2175 to 42% DR at Re = 10500. DR has been found to increase with A 
(we have measured up to 60% for A — X.2Up). The global energy budget has also been 
addressed: the energetic cost for producing the traveling waves is nearly minimum where 
DR is maximum, which implies large potential net savings. 

Backward-traveling waves always yield DR, while forward-traveling waves may also 
result in significant drag increase (DI) at finite phase speeds. Small levels of DI have 
also been observed at very low frequencies (infrared DI) for the case of uniform spanwise 
oscillations (n x = 0). 

Flow visualizations and turbulence statistics have been useful to show that, when 
DI takes place, the flow structures are tilted significantly along the spanwise direction 
and present a marked streamwise modulation, whereas they remain aligned with the 
streamwise direction when DR occurs. 

DI is observed when the waves travel with a phase speed comparable to the convection 
velocity of near-wall turbulent fluctuations. The near-wall structures are cyclically tilted 
by forcing, which ultimately results in an increased drag. At different phase speeds, DR 
is produced with a mechanism analogous to that of the oscillating wall. The optimal 
value of the oscillation period for the oscillating-wall regime still dictates the optimal 
conditions for the traveling waves, but the relevant period is that observed in a reference 
frame that convects with the flow. 
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